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We present a theoretical study of the spin transport properties of mono-atomic magnetic chains 
with a focus on the spectroscopical features of the I-V curve associated to spin-flip processes. Our 
calculations are based on the s-d model for magnetism with the electron transport treated at the 
level of the non-equilibrium Green's function formalism. Inelastic spin-flip scattering processes are 
introduced perturbatively via the flrst Born approximation and an expression for the associated self- 
energy is derived. The computational method is then applied to describe the I-V characteristics and 
its derivatives of one dimensional chains of Mn atoms and the results are then compared to available 
experimental data. We find a qualitative and quantitative agreement between the calculated and 
the experimental conductance spectra. Significantly we are able to describe the relative intensities 
of the spin excitation features in the I-V curve, by means of a careful analysis of the spin transition 
selection rules associated to the atomic chains. 



I. INTRODUCTION 

Inelastic electron tunneling spectroscopy (lETS) has 
become an important tool for investigating the elemen- 
tary excitations of nanoscale systema^. An excitation 
manifests itself as an abrupt change in the differential 
conductance of a two-probe device as the voltage sweeps 
across the excitation energy. This is the result of the 
opening of a new inelastic transport channel for the elec- 
tron tunneling. As the energy of the probed excitation 
sets the critical voltage and the temperature range where 
to observe lETS, it is not surprising that the first exper- 
iments were all related to molecular vibrations of rela- 
tively high energj0. 

More recently the continuous advances in low- 
temperature scanning tunneling microscopy (STM) have 
enabled the detection of excitations of different origin, 
namely those related to the spin degree of freedom. This 
different type of spectroscopy is usually named spin-flip 
lETS (SF-IETS). The first measurements of SF-IETS 
were for single atoms randomly deposited on surfaced. 
However, STM techniques also open the possibility of as- 
sembling and manipulating entire nano-structureP^ and 
of positioning magnetic ions on a surface at a desired 
position. This enables the construction of atomic mag- 
netic nano-structures and the study of the complex mag- 
netic excitations, resulting from the exchange interac- 
tion between the magnetic ions and the substrate, and 
also between the magnetic atoms themselves. Such a 
novel fabrication capability has produced a surge of ex- 
perimental studies on the spin excitations of magnetic 
nanostructures. These comprise the investigation of 
the conductance spectra of individual atomic spins^, of 
Kondo screened magnetic dimers^ and of multiple chains 
of Mn atom^. Interestingly similar investigations have 
now been extended to molecular chains composed of Co- 
phthalocyanines--. 

The recent rapid growth in the experimental activity 
has been matched by an equally fast explosion of theoret- 
ical works. A general and now standard approach to cal- 
culating the conductance spectra of the various possible 



magnetic nanostructures is that of combining a master 
equation solver for the quantum transport problem with 
model Hamiltonian describing the magnetic interactiorP. 
This is an intrinsic many-body approach, which in prin- 
ciple contains all the ingredients needed for solving the 
problem, once the various transfer rates are known. As 
such it usually requires a large number of parameters to 
be predictive. An alternative strategy consists in treat- 
ing the problem at the single particle level, by using a 
tunnel Hamiltonian and/or the standard Tersoff-Hamann 
descriptioiPSlfor the STM currenlfiiHll. This second class 
of computational scheme appears more amenable to be 
implemented together with first principles electron trans- 
port methods. However it still remains a hybrid theory, 
where the dynamic part of the problem needs to be ap- 
proached at the model Hamiltonian level, although the 
level of sophistication may include effects relate d to c ur- 
rent generated non-equilibrium spin population d^^ l ^^ l 

A standard theoretical approach to quantum transport 
is represented by the n on-equilibrium Green's function 
(NEGF) formalisnF^, whose mean-field version can be 
combined with state of the art electronic structure the- 
ories to produce efficient and predictive algorithm j^^ l ^° l 
This is to date the only fully quantitative quantum trans- 
port approach capable of scaling to large systemJ^U^ j e. 
it is the only one capable of performing simulations for 
materials set of current technological relevance. Impor- 
tantly for this discussion inelastic contributions to the 
elastic current can also be included within the NEGF for- 
malism. In the case of scattering to phonons the problem 
is usually treated perturbatively by constructing an ap- 
propriate self-energy at the level of either the first (IB A) 
or the self-consistent Born approximation (SCBA)^^. A 
similar approach to the case of spin excitations is cur- 
rently not available. 

There are two main reasons for this gap. On the one 
hand, the adiabatic separation between the electronic de- 
grees of freedom and those responsible for the inelastic 
excitations are well defined in the case of nuclear motion 
(phonons) but less clear in the case of spins, since even 
extremely localized spins have full electronic origir!^. On 
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the other hand, it is also unclear whether the perturba- 
tive approach is valid for spins, i.e. whether a suitable 
expansion parameter can be found. As such, as far as we 
know, an expression for a self-energy describing inelastic 
spin-flip events has not been derived so far. 

In the present work we undertake this challenge and 
formulate a perturbative theory of spin-flip spectroscopy 
based on the NEGF formalism for electron transport. 
Our theoretical analysis is based on a tight-binding 
Hamiltonian for the transport electrons, which are locally 
exchange coupled to quantum spins. This is essentially 
the s-d model for magnetisnl^^. We then proceed in con- 
structing an appropriate self-energy for the spin-degrees 
of freedom at the IBA level (note that for this particu- 
lar problem little differences are introduced by extending 
the treatment to the SCBA) and use this in the stan- 
dard NEGF scheme for transport. Our methodology is 
then applied to describing SF-IETS in atoms and atomic 
chains and results are favorably compared with experi- 
ments. 



II. INTERACTION HAMILTONIAN AND THE 
SCATTERING SELF-ENERGY 

The typical experimental setup considered in this work 
is that of an STM measurement, i.e. it comprises an 
STM tip positioned above one of the atoms of a mag- 
netic nanostructure, which in turn is weakly coupled to 
a metallic substrate across an insulating barrier. We 



made of the same material), which we describe by a one 
orbital per site tight-binding Hamiltonian with nearest 
neighbour interaction. 

The scattering region consists of a one-dimensional 
chain of N atoms. Each ith atom carries a quantum 
mechanical spin Si and it is characterized by an on- 
site energy Sq- We assume that the tip and substrate 
can only couple to one atom at the time in the scat- 
tering region and we describe such a coupling by the 
hopping integrals 7tip-s and 7sub-s, respectively. This 
means that the electronic states of the scattering region 
are broadened by the interaction with the electrode by 
Ttip-s/sub-s = 7tip-s/sub-s/7o where 70 is hopping pa- 
rameter in the leads. We assume that Eq ^ Ft/s lead- 
ing to a constant density of states at the Fermi energy. 
The electronic effect of the other atoms in the scattering 
region (i.e. creation of bonding and anti-bonding levels) 
can be neglected since these levels will also be far enough 
from the Fermi energy to ensure a constant density of 
states. 

The Hamiltonian of the scattering region contains 
three parts Hg = He + i?sp -I- i?o-sp, where Hq is the 
tight-binding electronic part, Hgp is the spin part and 
Hasp describes the electron-spin interaction. More ex- 
plicitly these three components write respectively as 
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FIG. 1: (Color online) Schematic representation of the device in- 
vestigated in this work. A scattering region, comprising spin- 
carrying atoms (light-blue circles) and described by the Hamilto- 
nian Hg , is sandwiched in between two semi-infinite electrodes (red 
circles). These mimic the substrate and the tip in a typical STM 
experiment. The electrodes are non-spin polarized and they are 
described by the Hamiltonian -ffguj, and //tip ■ In the scattering re- 
gion the transport electrons are exchange coupled to local quantum 
spins Si. 

model this system by using a pair of non-interacting semi- 
infinite leads that are separated by a scattering region as 
outlined in Fig. [T] The left-hand side lead, the scatter- 
ing region and the right-hand side lead represent respec- 
tively the STM tip, the spin-coupled nanostructure and 
the substrate and they are described by the Hamiltonian 
^tipj ^^s aiid -ffsub- For simplicity we assume an identi- 
cal electronic structure for both the leads (i.e. they are 
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Hc—sp ~ Jsd ^ ^ (c|^Q [""ilaa'Cia' ) ' Si . (3) 



The electronic part consists only in an on-site poten- 
tial, i.e. we neglect electron hopping between the sites 
[the electron ladder operators cl^/cia create/annihilate 
an electron at site i with spin a (=t, J,)]. 

We model the spin-spin interaction between the local- 
ized spins {Si} by a nearest neighbour Heisenberg Hamil- 
tonian with coupling strength Jm- Furthermore we in- 
clude interaction with an external magnetic field B (/ie 
is the Bohr magneton and g the gyromagnetic ratio) and 
both uni-axial and t ransv erse anisotropy of magnitude 
D and E respectivelj^^H xhe electron-spin interaction 
Hamiltonian is constructed within the s-d modeP^ where 
the transport electron, s, are locally exchanged coupled 
to quantum spins, {S;} {d indicates that the local mo- 
ments originate from the atomic d shell). In this case 
the interaction strength is Jsd and in equation ([s]) cr is a 
vector of Pauli matrices. 

We now proceed to deriving the perturbation theory in 
the electron-spin interaction (note that He and Hgp in- 
dependently can be easily diagonalized) . Our strategy is 
that of first constructing the many-body Green's function 
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at the level of IBAt^SHSU and then, by Dyson's equation, to 
evaluate the interacting self-energy. In particular we fol- 
low closely the procedure laid out in reference [53] ■ Since 
we consider a non-equilibrium situation at zero tempera- 
ture our starting point is the Keldysh^^ contour-ordered 
single-body Green's function in the many-body ground 
state. 



[G{Ty)U, = -i{\Tc{c„{r)cUr')}\), (4) 



where the time-average is performed over the full in- 
teracting ground state |) (note that for clarity we have 
dropped the site index z, which will be explicitly included 
only where necessary). Equation Q can be expanded up 
to the n-th order in the interaction Hamiltonian, Hcsp, 
as 
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(0|rc{ge-sp(ri) . . .ge-sp(r„)cfc^(r)4,,,(rO}|0) 
U (— oo, — oo) 



where U is the time-evolution unitary operator and the 
time-averages are now over the known non-interacting 
(Jsd — 0) ground state |0). The time integration over 
T is ordered on the contour C going from — oo to -l-oo 
and then returning from +oo to ~oo, since the ground 
state of the non-equilibrium system can only be defined 
at -ocP. If the expansion is truncated to the first order 
one obtains a Zeeman-like term which can be neglected 
as long as Eq ^ 7tip-s/sub-s- The first contribution of 
interest then appears at the second order. This can be 
obtained by inserting the explicit expression for -ffc-sp {t) 
[Eq. ^] into Eq. Q 



(0|rcK(r)cl(Ti)c„,(Ti)4(r2)c^.(T2)c^,(r')}|0) 
X Wc{S\r^)S^{T:,)m\A.A'y%!i'. (6) 

where the indices i and j now run over the cartesian 
coordinates a;, y and z. 

A full contour-ordered expansion must now be per- 
formed on both the electron bracket and the spin bracket. 
The electron bracket has six different time-ordering com- 
binations, which are explicitly listed in the equation y\ 



below 

(0|TcK(r)cl(Ti)c„,(Ti)4(r2)c^.(T2)4,(r')}|0) 
= (0|rc{c.(T)4(ri)}|0) X (0|Tc{c„,(ri)4(r2)}|0) 

x(O|Tc{c0,(T2)4,(r')}|O) 
+ (0|TcK(r)4(T2)}|0) X (0|Tc{c„,(Ti)4,(r')}|0) 

x(0|Tc{c„(T2)4,(ri)}|0) 
+ (0|rc{c.(T)4(ri)}|0) X (0|Tc{c„,(Ti)4,(r')}|0) 

x(0|Tc{c^(T2)c^,(r2)}|0) 
+ (0|Tc{c,(r)4(T2)}|0) X (0|Tc{c^j,(T2)4,(r')}|0) 

x(0|rc{cl(ri)c„-(Ti)}|0) 
-H(0|rc{c.(T)4,(T')}|0) X (0|rc{cl(ri)c„,(ri)}|0) 

x(0|Tc{4(T2)c^,(r2)}|0) 
-(0|Tc{c.(r)4,(r')}|0) x (0|Tc{c„, (ri)4(r2)}|0) 

x(O|rc{c0'(T2)ct(ri)}|O). (7) 

The last two combinations can be eliminated since they 
represent unconnected Feynman diagrams which vanish 
in the averaging process^Sl. The first and the second term 
are equal under the exchange of the indexes, and so are 
the third and the fourth ones. This leaves us with a 
simplified expression which, when compared to Eq. ([5|, 
gives us 

(0|Tc{c.(T)cUri)c„,(n)4(r2)c^'(^2)4,(T')}|0) 

= 2i^5aa5a'P^li'a'\GQ{T,Txy\aa\GQ{Tx,Ti)\aia' 
X {Gq{t2,t'%,„, 

-I- It'^aa^a'a'^fili' [Go(t, Ti )] [Gq (^1 , t' )] ct'ct' 

X [Go(t2,T2)]/3;3. (8) 

In this case, since the averaging bracket is over the non- 
interacting ground state. Go represents the un-perturbed 
electronic Green's function and can be calculated exactly. 

We then return to Eq. ([6]) and evaluate the spin 
bracket. The ground state of the non-interacting spin 
system alone can be found by diagonalizing exactly -ffgp- 
This is achieved by constructing the full basis { | n) } where 
n = —S, —S + 1, +S. The resulting eigenvectors, |m), 
and eigenvalues, E„i, satisfy the Schrodinger equation 
Hsp\m) — em|m) and they can be used to re-write the 
operators S''(r) as 

S\r)=Y,{m\S^\n)dl{r)d„{T). (9) 

m,n 

Here c?„ is a desctruction operator for a quasi-particle of 
the spin system. The quasi-particles are then assumed to 
be fermionic in nature so that they obey the anticommu- 
tation rules [d}„,c?„] = 6mn and [d3„:4J = [dm,dn] = 0. 
Such an assumption is valid as long as the excitations 
considered are always around the ground state, i.e. under 
the condition that the spin system can always efficiently 



4 



relax back to the ground state between spin-flip events. 
We can then define a contour-ordered spin Green's func- 
tion as follows 



[D{t,t')] 



^^{\Tc{d^{T)dl{r')}\) . (10) 



By inserting the expressions in the equations (|9| and ( 10 ) 
into the spin bracket and by computing the time-ordered 
contraction we finally obtain 

(0|Tc{5+(ri)57(r2)}|0> = 

(11) 



The set of equations \m and ( 11 ) can now be incorpo- 



rated into the expression for the second order contribu- 
tion to the many-body Green's function [Eq. ([6|]. Then, 
by using Dyson's equation^^, one can finally write the 
second order contribution to the interacting self-energy 
which reads 



X [Go(ri,r2)]^^EH^»("l^'k> 

(12) 

If we now assume that the ground state electronic spin 
levels are degenerate, i.e. [GoJ-f^ = [Go]^^, then the 
only quantity of interest is the trace of the self-energy 
over spin indices. By performing such a trace, the spin- 
independent self-energy finally reads 

Si„t(ri,r2)(2) = ^ |(™|5>)pGo(ri,r2) 

(13) 



where we have used the results Tr[(TV-'] = 5ij and 
Tr[cr*] = 0. 

At this point we can calculate the real-time quantities, 
such as the lesser (greater) self-energies, by using the 
Langreth's theorem for time ordering over the contour, 
Ti e Gi(G2) and T2 G C2(Gi)^0'. Gi is the time ordering 
contour from —oo to +oo and G2 is the time anti-ordering 
contour from -foo to —00. We find 

^Ut„t,)^^^=-2Ji^ J2 |(m|5»pGf(ti,t2) 

i.m.n 

X [D'^{ti,t2)]n..n[DQ{t2,ti)] 

= -2j2, ^ \{m\S^\n)\^Gt{hM) 



X P„(l-F™)e 



z,m,n 
±i(e™-e„)(ti-t2) 



(14) 



where in the second step we have written explicitly 
D^{ti,t2) in terms of the spin-level occupations, P„ = 



d\^dn- The dependence of I]?^ over the energy, E, can 
found by simple Fourier transform 



(2) 



-2JI^ E \{MS»?Pn{l~Prn) 



X G^{E±Vt,, 



(15) 



where 0.mn = £m — £n and where the + (— ) symbol cor- 
responds to (^^)- 

Electron-spin scattering even ts a re now fully described 
by [Efj^^{E)]. In particular Eq. (15) describes the process 



where an incoming electron with energy E experiences 
a spin- flip process, which changes its energy by ±f2m„. 
This is the result of the electron-spin interaction with the 
local spins. Such a process is schematically represented 
in Fig. [2j Note that the probability for an excitation to 
occur is determined by the prefactors {m\S'^\n)\'^Pn{l — 
Pm), i-e. by the state of the spin system. 




incoming e 



spin system 



outgoing e 



FIG. 2; Schematic representation of the inelastic process described 
by the greater self-energy, S?^^. An incoming electron scatters 
against a localized spin and decreases its energy by Hmn- This is 
transfer to the local spin system, which undergos a spin transition 
|n) ^ \m). 
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III. NON-EQUILIBRIUM GREEN'S FUNCTION 
METHOD FOR ELECTRON TRANSPORT 

The transport properties of the device are described 
by using the non-equilibrium Green's function (NEGF) 
method. This has been extensively described in the 



past! ^'"^" l Here we only summarize the main concepts 
and we highlight the modifications associated to describ- 
ing electron-spin interaction. A two-probe device can 
be divided into three distinct regions, two semi-infinite 
leads which represent the STM tip and the underlying 
substrate and a central scattering region (see Fig. [T]). 
The leads act as charge reservoirs and they are charac- 
terized by their chemical potential, respectively /xtip and 
/isub- The external bias is introduced in the form of a 
relative shift (symmetric) of the two chemical potentials 
with respect to each other. The underline assumption 
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of the method is that under the external bias there is 
no rearrangement of the electronic structure of the leads, 
i.e. that electron screening in the leads is efficient. This 
simplifies the problem to that of calculating the retarded 
Green's function of the scattering region^^ only 



G{E) ^ lim [{E - iri)I - - ^{E)Y 



(16) 



Here He is the electronic part of the Hamiltonian [Eq. ([ij] 
and is the retarded self-energy, which incorporates 

the effects of the leads and of the inelastic interaction. In 
particular 'S(E) writes as 



(17) 



where Stip(-E) and T,suh{E) are respectively the STM tip 
and substrate self-energies, while Sjnt is the scattering 
self-energy describing the electron-spin interaction. The 
leads' self-energies can be written in terms of the surface 
Green's functions ^tip and ^sub and the coupling matri- 
ces between the leads and the scattering region Has 
(a=tip, sub) 



tip 



-s 9tip{E) i?tip-s 



Ssub(£^) — -f^sub-S 9suh{E) -f^J^b-S 



(18) 
(19) 



The surface Green's functions can be found by first 
constructing the Green's function for an infinite sys- 
tem and then by applying the appropriate boundary 
conditions'^. For a single-site nearest-neighbour one- 
dimensional tight-binding chain ga takes a simple close 
form 



1 



gtipiE) = gsuhiE) = — e 
70 

E~eo 



ik(E) 



k{E) = cos~^ 



270 



(20) 
(21) 



where we k[E) and 70 are respectively the dispersion 
relation and the hopping parameter. We assume equal 
hopping parameter in the leads and also the condition 
7o > £0- The interaction between the scattering region 
and the leads has the effect of broadening the scattering 
region discrete energy levels. The broadening functions 
write 



T„{E)^i[^„{E)-Y.l{E) 



(22) 



where (a=tip, sub). The resulting lesser/greater leads 
self-energies are related to the broadening and to the 
population of the electrons/holes in the leads 



1:<{E)^ f^{E)T^{E), 
Y>{E) = [l- UE)]T^{E), 



(23) 
(24) 



where ^(S) = 1/ 
tions. 
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cxp 



are the Fermi func- 



Finally, the retarded electron-spin scattering self- 
energy is found from the Hilbert transfer 



dE' i:UE')-^UE') 
27r E' -E 

,^UE)-^UE) 



(25) 



where W denotes the principal value and where the full 
expressions for Yif^^i^E) and Yi'^^^^{E) have been given in 
the previous section. Note that at variance with a similar 
expression for the self-energy describing electron-phonon 
scattering, in equation ( 25 ) there is no Hartree-like con- 



tribution. This usually describes polaronic effects and it 
is known to cause the breakdown of the Born approxi- 
mation for strong electron-phonon coupling'^- . 

Additional quantities of interest, which enter in the 
definition of the electron current, are the lesser and 
greater Green's functions (G< and G>). Their diagonal 
elements respectively give the electron and hole density 
in the scattering region. These are defined by 



G^(£;) = G{E)[Y.'>{E)]G^{E) , 



where 



Y.^{E) 



^UE) 



(26) 



(27) 



As already mention the application of a potential dif- 
ference V across the device produces a shift in the tip 
and substrate chemical potentials. By assuming, as done 
here, that the tip and the substrate share the same 
Fermi level Ep then we have /^tip ^ Ep + eV/2 and 
A'sub = Ep + eV/2. This is equivalent to replacing the 
energy in the leads with E — )■ £'±6^/2 (e is the electron 
charge) . Finally the current can then be calculated for a 
range of values of V at any of the leads a as 



+00 



Ia{E) dE , 



(28) 



' —00 



UE) = -Tr{[i:>{E)G<{E)] - [£< (i?)G> (i?)]} . (29) 



IV. RESULTS 

The entire procedure outlined in the previous two sec- 
tions is now implemented numerically in order to calcu- 
late the I-V characteristics and the conductance spectra 
{dl /dV-V) of some selected magnetic nanostructures. In 
particular we aim at reproducing the dl /dV-V curve ob- 
served through Mn monoatomic chains deposited on thin 
CuN insulating atomic layers, as described in reference 
[7] . Clearly our ID model does not include all the details 
involved in the actual SF-IETS STM experiments, but 
we argue here that it contains already all the ingredients 
to reproduce the main experimental features. 
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Quantity 


Symbol 


Value 


Origin 


[Ref 


Atomic Apin 


s 


5 
2 


Exp. 


111 


TcmpsrcLturs 


T 


0.6 K 


Exp. 


Ui 


d-d exchcinge 


J dd 


+6.2 meV 


Exp. 


[71 


s-d excliange 




+500 meV 


DFT 


l33l 


Fermi Energy 


£p 


meV 


DFT 


l33l 


Lead hopping integral 


'Yn 


10000 meV 


DFT 


031 


Channel on-site energy 


£o 


1000 meV 


DFT 


B31 


Tip-channel hopping 


7tip-S 


50 meV 


DFT 


m 


Substrate-channel hopping 


7sub-S 


500 meV 


DFT 


m 


Axial anisotropy 


D 


-0.037 meV 


Exp. 


u 


Transverse anisotropy 


E 


0.007 meV 


Exp. 


m 



TABLE I; Empirical parameters used in the numerical simulations 
presented in this work and their assigned numerical values. In the 
rightmost column we report the source of a particular estimate. 
"Channel" here means the scattering region. 

In Table I we list all the parameters used in our sim- 
ulations and their assigned values. These have been ei- 
ther inferred from the experiments (Exp^ or have been 
estimated from density functional theory (DFT) calcula- 
tion^. The local Mn spin is set to be 5/2, as proposed 
in the original experimental work^^, confirmed by DFT 
calculationj22l and expected from the nominal Mn va- 
lence. The spin-spin coupling parameter Jdd corresponds 
to an antiferromagnetic order between the neighbouring 
Mn spins, a feature verified in the experimental conduc- 
tance spectra. The lead on-site energy is suitably set to 
zero and simply defines the reference potential. We also 
note that the value for Jsd is determined from theory 
to be of the order of SOOmeV'^S. Although we employ a 
zero-temperature formalism we nevertheless evaluate the 
Fermi functions of the leads at the small temperature of 

0. 6 K. This allows us to include minor thermal smearing 
of the electron gas in the leads and consequently of the 
conductance profile. Finally we notice that the scatter- 
ing region is expected to be significantly more strongly 
coupled to the substrate than to the tip. This essentially 
means that the spin population in the scattering region 
remains effectively in equilibriuirP^ with the substrate, 

1. e. there is no accumulation of spin excitations between 
two inelastic tunneling events. 

Figure [3] shows the calculated conductance spectra 
(normalised against the elastic contribution to the con- 
ductance Go) for A''-atom long Mn chains [N < 4) in no 
external magnetic field. In the case of iV=:3 and N—A the 
spectrum is calculated for the STM tip placed above the 
second atom in the chain. From the figure it is clear that 
all the qualitative experimental findings are correctly re- 
produced. In particular we notice the strong dependence 
of the conductance profile over the parity of the chains, 
with chains comprising an odd number of atoms (odd 
chains) exhibiting a conductance dip at around ^=0, 
which is absent for even chains. 

The ground state of each odd chain has a net total spin 
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FIG. 3: (Color online) Normalised conductance spectra for Mn 
chains of different lengths A'^. The various spectra are offset for 
clarity. The tip is placed above the second atom of the chain for 
chains with N > 2. We notice the strong dependence of the spectra 
on the chain parity. 

5'tot = 5/2. This is affected by the transverse and axial 
anisotropy, which lift the ground state degeneracy and 
allow a transition between the ground state and the first 
spin-flip excited state to occur. The excitation of such a 
transition results in a conductance step at a voltage cor- 
responding to the transition energy. Since the anisotropy 
energies are small (see Table I), the excitation energy is 
small as well and the feature in the conductance profile 
appears near V—0. In contrast the even chains do not 
show any zero-bias excitation, since the ground state is 
a singlet and does not carry a magnetic moment. This 
is a direct proof that the magnetic interaction between 
Mn ions is antiferromagnetic. Should this have been fer- 
romagnetic, even chains would also have shown zero-bias 
anomalies. 

As mentioned above the conductance steps encoun- 
tered at around V=0 correspond to spin-flip events, i.e. 
to electron scattering processes that produce the transi- 
tion |m = 5/2) — > |m = 3/2) (m is the magnetic quantum 
number) but that also preserve the total spin ^tot = 5/2. 
The first net spin changing transition occurs for N — 2 
and corresponds to the large conductance step found at 
V — Jdd/e = 6.2 mV. This is investigated in more de- 
tail in figure [4j where we also include the dimer's spec- 
trum calculated when a magnetic field is applied along 
the z-direction. The effect of the magnetic field is that of 
splitting the single excitation line into three distinct con- 
ductance steps. In this case in fact the antiferromagnetic 
dimer has a singlet {Stot=0) ground state and a triplet 
(5'tot=l) first excited state. Therefore, an excitation from 
the ground state to the first excited state changes the net 
spin of the dimer. When a magnetic field is applied the 
degeneracy of the triplet excited state is lifted and exci- 
tations having the three Zeeman-splitted levels as final 
states are possible. This produces the three-fold split- 
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ting of the conductance steps as seen to agree well with 
experiment in figure |4] 
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FIG. 4: (Color online) Theoretical (A) and experimental (B) con- 
ductance spectra for the N=2 chains in a finite external magnetic 
field. For B ^ we note a splitting of the conductance steps, cor- 
responding to the Zeeman split of the final triplet excited state. 
This is observed in the experiments of reference [7], which are also 
reported in panel (B). 



Figure [5] shows the conductance spectrum for the 
trimer (A'^=3). This chain exhibits similar trends as 
those found for the dimer as a magnetic field is applied, 
i.e. there is a Zeeman split of the zero-field conductance 
steps. From the figure one may notice that for both the 
theoretical (A) and experimental (B) data there is a shift 
of the principle step to lower voltages as the magnetic 
field is increased. We remark here that in order to recre- 
ate the conductance profile of the trimer, additional ferro- 
magnetic second-nearest-neigiibour interaction between 
the local spins at the edge of the chain must be included 
m the modeP. The magnitude of this additional ex- 
change parameter is approximately half of that of Jdd- 
The inclusion of such a second-nearest-neighbour cou- 



pling constant changes the position of the conductance 
step from a second excited state at 27 mV to a first ex- 
cited state at 16 mV. This correction is also included in 
the calculations for 7V=4 (see Fig. |3|, again giving good 
agreement with experiments. 

In all of the spectra investigated the most striking 
agreement with the experiments concerns the correct pre- 
diction of the SF-IETS excitation voltages. In particu- 
lar not all the possible spin excitations feature in the 
AI /dV — V curve as a result of spin-selection rules. A 
careful analysis of the inelastic process, as outlined in 
the first section, reveals why some excitations occur and 
why some other do not. In particular we note that the 
proper selection rules for transitions are gove rned by the 



prefactor in the self-energy |(m|S'*|n)p [Eq. (15)]. This 



effectively selects which excitations are more probable 
to occur. For example, the full energy spectrum of the 
trimer has 6^ = 216 eigenvalues but only a small por- 
tion of these eigenvalues will contribute significantly to 
the scattering self-energy. Consequently only a few tran- 
sitions will have influence on the conductance spectrum. 
This for instance explains why the first conductance step 
of the dimer is considerably larger than that of other 
chains of different lengths. In fact our calculations show 
that a transition from the singlet ground state {S=Q) to 
any of the triplet excited states (S'=l; m= -1, 0, -1-1) 
has equal probability. This results in a conductance step 
which is approximately three times larger than any single 
spin-flip event that occurs in an odd chain. 



V. CONCLUSIONS 

In conclusion we have formulated a theory of elec- 
tron transport through magnetic nanostructures based 
on the NEGF formalism and including inelastic spin- 
scattering. We have first explicitly derived an expres- 
sion for the interacting self-energy at the level of the 
first Born approximation in the electron-spin interaction. 
This has been used to calculate the current-voltage and 
the conductance-voltage curves for a ID system of mag- 
netically coupled spins at finite bias and in a magnetic 
field. Our results reproduce well the features of recent 
SF-IETS experiments for ID Mn chains. Most notably 
the severe dependence of the conductance spectra on the 
chain parity and the selection-rule suppression of certain 
transitions is a direct outcome of the theory. Our NEGF 
approach is therefore a valid alternative to master equa- 
tion based schemes, with the advantage that it can be 
scaled up to larger systems and combined with more so- 
phisticated electronic structure theories. 
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FIG. 5: (Color online) Theoretical (A) and experimental (B) eonductance spectra for the N = 3 trimer at different magnetic field. We 
notice the shift of the first conductance step transition to lower voltages as the magnetic field increases. This is in good agreement with 
STM experiments of reference 0, which are reproduced here in panel (B). 
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